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ABSTRACT 

The accretion-powered millisecond pulsar ICR J00291+5934 underwent two ~10 d long outbursts 
during 2008, separated by 30 d in quiescence. Such a short quiescent period between outbursts has 
never been seen before from a neutron star X-ray transient. X-ray pulsations at the 599 Hz spin 
frequency are detected throughout both outbursts. For the first time, we derive a pulse phase model 
that connects two outbursts, providing a long baseline for spin frequency measurement. Comparison 
with the frequency measured during the 2004 outburst of this source gives a spin-down during qui- 
escence of —(4 ± 1) x 10~ 15 Hz s" 1 , approximately an order of magnitude larger than the long-term 
spin-down observed in the 401 Hz accretion-powered pulsar SAX J1808.4— 3658. If this spin-down 
is due to magnetic dipole radiation, it requires a 2 x 10 8 G field strength, and its high spin-down 
luminosity may be detectable with the Fermi Large Area Telescope. Alternatively, this large spin- 
down could be produced by gravitational wave emission from a fractional mass quadrupole moment 
of Q/I = 1 x 10~ 9 . The rapid succession of the outbursts also provides a unique test of models 
for accretion in low-mass X-ray binaries. Disk instability models generally predict that an outburst 
will leave the accretion disk too depleted to fuel a second outburst after such a brief quiescence. We 
suggest a modification in which the outburst is shut off by the onset of a propeller effect before the 
disk is depleted. This model can explain the short quiescence and the unusually slow rise of the light 
curve of the second 2008 outburst. 

Subject headings: binaries: general — stars: individual (IGR J00291+5934) — stars: neutron — stars: 
rotation — X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

The longevity of the Rossi X-ray Timing Explorer 
(RXTE) has provided the opportunity to observe multi- 
ple outbursts from accretion-powered millisecond pulsars 
(AMSPs) with recurrence times of ~10 yr or less. Such 
observations can address a diverse array of science. The 
pulse timing of multiple outbursts can measure changes 
in the spin of the neutron star (NS), which places lim- 
its on its magnetic field and gravitational wave emission. 
Comparison of the light curves and spin frequency deriva- 
tives during outburst probes the interaction between the 
NS and the accretion disk, while the shape of the pulses 
constrains the nature of the X-ray emission and the mag- 
netically channeled accretion flow. Finally, the observa- 
tion of multiple outbursts provides tests of disk insta- 
bility models and other theories put forward to explain 
the recurrence of these sources. In this paper, we report 
on the observation by RXTE of a second and third out- 
burst of the 599 Hz accreting pulsar IGR J00291+5934. 
These outbursts were separated by only 30 d of quies- 
cence, a more rapid recurrence than has ever been seen 
before from a NS low-mass X-ray binary (LMXB). The 
proximity of the two outbursts proves particularly useful 
for addressing many of the above questions. 

The International Gamma-ray Astrophysics Labra- 
tory {INTEGRAL) first detected IGR J00291+5934 
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at the onset of a n outburst on 2004 December 2 
(jEckert et al.l 12004ft . Follow-up observations with 
RXTE revealed pulsatio ns at a frequency of 598.89 Hz 
(iMarkwardt et al.ll2004bf) . mod ulated by a 147.4 minute 
orbit (Markwa rdt et al.ll2004aft . Analysis of the pulsa- 
tions revealed them to be highly sinusoidal across a wide 
range of energies ( Ga lloway et al.l 120051 : iFalanga et al.1 
l2005bft . Pulses arrived progressively sooner with increas- 
ing energy over 2-8 keV, followi ng the pattern of sof t 
lags observed in other AMSPs (jGallowav et al.l l2005j) . 
but above 8 keV these soft lags diminished, a reversal not 
seen in other AMSPs that may ha ve important implica- 
tions for the origin of these lags (jFalanga fc Titarchukl 
12007ft . Fractional amplitudes were be tween 5%-10% 
rms, generally decreasing with energy (j Galloway et al.l 
2005). Pulse timing models for the 2004 outburst 
require a spin deriva t ive of 8.5(1.1 ) x 10 ~ 13 Hz s _1 
(jFalanga et all l2005bl: iBurderi et all 12007ft . which has 
been ascrib ed to the NS being spun up by the accret- 
ing matter (|Burderi et a l. 2007). (Note that parentheti- 
cal uncertaintes are at the 1 a level throughout this pa- 
per.) Aperiodic timing of IGR J00291+5934 reveals an 
unusual amount of timing noise at very low frequencies 
(0.01-0.1 Hz), producing a timing spectrum more akin 
to black holes than to other NS low-mass X-ray binaries 
(|Linares et al.ll2007ft . 

A mission-long light curve from the RXTE All Sky 
Monitor showed marginal (5cr) detections of earlier 
outbursts during 1998 November and 2001 September 
(jRemillardl 12004ft . With each outburst, the duration 
of quiescence increased by 160-170 d. A quadratic 
fit to these outburst times predicted a 3.6 yr quies- 
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cence between t he 2004 outbur st and the start of the 
2008 outbursts (|GalIowavl 120081) . This estimate proved 
to be accurate to within 1% of the recurrence period, 
with the source returnin g to outburst on 2008 Aug 13 
(jChakrabartv et al.ll2008h . 

X-ray observations between the 2008 outbursts show 
that IGR J00291+5934 reached a flux level consistent 
with prior measurements during longer intervals of qui- 
escence. On Aug 21 the Swift X-ray Telescope (XRT) 
gave a 3 a upper limit on the unabsorbed 2-10 keV flux 
of 4.7 x 10~ 12 erg cm~ 2 s _1 , and on Aug 25 an XMM ob- 
servation detected the source at an unabsorb ed 2-10 keV 
flux of (1.4 ± 0.3) x 10" 14 erg cm" 2 s" 1 (|Lewis et al.l 
1201(1 . In comparison, a Chandra observation nearly two 
years after the 2004 outburst showed an un absorbed 0.5- 
10 ke V flux of 1 x 10~ 13 erg cm~ 2 s" 1 (jJonker et al.1 
2008). Because the accretion episodes during 2008 were 
separated by a quiescent period during which accretion 
was very low or halted entirely, we refer to them as two 
separate outbursts, 2008a and 2008b. 

In this paper, we present a detailed analysis of the 2008 
outbursts of IGR J00291+5934 and the implications for 
our understanding of recurrent X-ray transients. We de- 
scribe the RXTE observations and data analysis in §2, 
and we present the results of this analysis in §3. In §4 we 
examine what the 2008 double outburst can tell us about 
the theories of accretion disk instability. In §5 we con- 
sider other evidence for changes in the accretion regime 
at low fluxes, and in §6 we derive limits on the NS mag- 
netic field. In §7 we discuss the long-term spin evolution 
and the possible sources of torque on the NS. In a final 
section, we summarize our results and conclusions. 

2. OBSERVATIONS AND ANALYSIS 

IGR J00291+5934 was found to be in outburst dur- 
ing an RXTE monitorin g observation on 2008 Aug 13 
(jChakrabartv et al.1 l200l . RXTE observed the source 
daily during this outburst (observation IDs 93013-07-*). 
By Aug 22, the source had faded below the RXTE de- 
tection threshold (~1 x 10~ n erg cm~ 2 s _1 ). Intensive 
RXTE observation ended on Aug 30, and the program 
of twice- weekly «1 ks monitoring observations (obser- 
vation IDs 93435-01-*) resumed. On Sep 18, an opti- 
cal observation found the source to have unexpectedly 
re-brightened, and a Swift XRT observation on Sep 20 
confirmed that IG R J00291+5934 was again in outburst 
(|Lewis et a l. 2008). Longer, more frequent RXTE obser- 
vations resumed shortly thereafter. This second outburst 
persisted through Oct 3. During the first outburst, 20 
RXTE observations were taken, totaling 75.3 ks of ex- 
posure; during the second, 12 observations were taken, 
totaling 37.9 ks. Figure [T] shows the RXTE light curve 
of these outbursts and the times of observations. 

We analyzed the da ta from the RXTE Proportional 
Counter Array fPCA; fJahoda et al.l 119961 ). which com- 
prises five identical, co-aligned proportional counter units 
(PCUs) sensitive to 2.5-60 keV photons within «1° of the 
pointing axis. Due to the increase d frequency of high - 
voltage breakdowns in the PCUs (|Jahoda et al.l [2006). 
an average of only 1.7 PCUs were taking usable data 
during the observations. (For comparison, an average 
of 2.8 PCUs were active during the 2004 outburst of 
IGR J00291+5934.) The resulting mean effective area 
during the 2008 outbursts was 2100 cm 2 . 
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Fig. 1.— The RXTE PCA light curve of the 2008 outbursts from 
IGR J00291+5934. The horizontal dotted line shows the estimated 
contribution from the nearby intermediate polar V709 Cas. 

We derived outburst light curves using only PCU 2, 
which was active throughout all the observations. We fit 
each observation with an absorbed blackbody plus power 
law model using XSPEC and the latest PCA response 
matrix (version 11.7) to estimate the 2.5-25 keV flux 
and its uncertainty. Note that the absorbed fluxes are 
given throughout this paper (i.e., the flux as observed by 
RXTE). Instrumental background levels were estimated 
using the FTOOL pcabackestQ 

The large field of view of the PCA also admits photons 
from the nearby intermediate polar V709 Cas, located 17 
arcmin from IGR J00291+5934. For all observations the 
RXTE science axis was pointed at IGR J00291+5934, 
resulting in the PCA collimator admitting 70% of the 
photons from V70 9 Cas (using the l inear collimator 
model from §8.3 of Uahoda et al.ll2006l ). We estimated 
its contribution using the monitoring observations of 
the IGR J00291+5934 field that preceded and followed 
the 2008 outbursts. The mean 2.5-25 keV flux was 
0.68 x 10 -10 erg cm~ 2 s _1 after subtracting the instru- 
mental background, implying an unattenuated flux from 
V709 Cas of 1.0 x 10~ 10 erg cm" 2 s" 1 . Earlier obser- 
vations suggest that its 2.5-25 keV flux varies within 
(0.9-1.3) x TO- 10 erg cm" 2 s" 1 (|de Martino et all 12001 
iFalanga et al.ll2005cJ : and references therein), consistent 
with our results. The uncertainty in our V709 Cas 
flux measurements sets a PCA detection threshold of 
-1 x 10" 11 erg cm" 2 s" 1 for IGR J00291+5934. The 
flux of V709 Cas is modulated by its 312.8 s rotational 
period (jHaberl fc Motch1ll995D with a fractional ampli- 
tude of 20^SO%^Mm^m^MMoM)- When dividing 
RXTE observations into shorter intervals, we chose in- 
tegral multiples of the V709 Cas period to simplify the 
estimation of its contribution. 

For our coherent timing analysis, we included 2-15 keV 
photons to maximize the signal to noise ratio. We shifted 
the photon arrival times to the solar system baryccnter 
using the DE405 solar system ephe meris and the optica l 
and near-infrared position given bv lTorres et~aT1 ([2008): 
R.A. = 00 h 29 m 03 i !05±0 ! ;01, decl. = +59 o 34 / 18 / .'93±0 / . / 05 

5 http : //heasarc .gsf c .nasa. gov/docs/xte/recipes/pcabackest .html 
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Fig. 2. — The light curves, phases, and fractional amplitudes of the observed outbursts of IGR J00291+5934. The light curves show 
the 2.5—25 keV flux, with one point per observation, after subtracting the contribution from instrumental background and the nearby 
intermediate polar V709 Cas. The dotted gray line indicates the critical flux f x ,c, which marks the quickening of the decay during the 2004 
and 2008a outbursts and the maximum during the 2008b outburst. The pulse phases show the residuals measured relative to the best-fit 
2008 frequency, uo = 598.892130804 Hz. The gray lines show the best-fit models for each outburst. Phases and fractional amplitudes are 
shown for the 2.5-15 keV band. 95% upper limits on the fractional amplitudes are indicated for observations without detectable pulsations; 
these upper limits for the first and last observations of the 2008b outburst were very weak (~50% rms). 



(J2000.0). This position is consist ent with the Chandra 
X-ray position (iPaizis et al.l I2005D an d the earlier opti- 
cal position of iFox fc Kulkarnil (|2004| ) . which was used 
to der ive the pulse timing ephemeris in IGallowav et al.l 
(2005)0 After barycentering the photon arrival times, 
we applied the RXTE fine clock correction and filtered 
out data during Earth occultations and intervals of un- 
stable pointing. We searched the data for thermonuclear 
(type I) X-ray bursts; none were present. 

To measure the times of arrival (TOAs) and fractional 
amplitudes of the persistent pulsations, we folded 626 s 
intervals of data (twi ce the V709 Cas period) using the 
IGallowav et al.l (|2005l ) timing model. The timing models 
were applied and fitted using the TEMPO pulse timing 
program, version f f .005Q and assumed a circular orbit 
and a fixed spin frequency. From the resulting folded 
pulse profiles, we measured the phases and amplitudes 
of the fundamental harmonic of the pulsation. Higher 
harmonics had insufficient power to be useful for pulse 
phase timing. Estimation of phase and fractional am- 
plitude uncerta i nties follow the procedures described in 
lHartman et~aTl (|2008l ). 

We also analyzed the aperiodic variability during the 

6 This R.A. d i ffers b y -Of 032 (3.2 a) relative to the radio position 
of IRupen et al. (2004), which was us ed to derive the pu l se tim ing 
ephemerides of Falanga et al. (2005b) and Burdcri et al. (2007). If 
these ephemerides had used the optical position, their frequencies 
would shift by Ai/ = 4 X 10" 8 Hz and Au = 1.0 X 10~ 14 Hz s -1 . 
This offset in v is far too small to account for the spin-ups reported 
by the se authors. 

7 At http://wwww.atnf.csiro.au/rcsearch/pulsar/tempo 



IGR J00291+5934 outbursts to characterize its broad- 
band noise properties and search for quasi-periodic os- 
cillations ( QPOs). We fo l lowed the procedures and con- 
ventions of iLinares et al.1 (|2007t ). who measured the tim- 
ing spectrum of the 2004 outburst, so that our results 
would be directly comparable. We selected 2.5-30 keV 
photons, as per the earlier study, and modified their ar- 
rival times using our orbital ephemeris to shift them into 
the frame of the NS. We then performed Fourier trans- 
forms on 1024 s segments of data and normalized the 
resulting power sp ectra using the rms normalization of 
Ivan der Klisl ([1995). We estimated the V709 Cas contri- 
bution as previously described and corrected for it and 
the instrumental background. We averaged the power 
spectra to create a single spectrum for each outburst, 
then fit them with multipl e Lore ntzians following the 
conventions of iBelloni et al.l (|2002f ) . 

3. RESULTS 

3.1. Outburst Light Curves and Fluences 

The light curve of the first 2008 outburst showed a 
fast rise and slow decay. The rise of the outburst oc- 
curred in less than 3.5 days, the length of time sep- 
arating the observations during the PCA monitoring 
campaign. A monitoring observation on MJD 54688.4 
gave a 2.5-25 keV flux from the field of (0.72 ± 0.06) x 
10" 10 erg cm" 2 s _1 , consistent with the flux of V709 Cas. 
By the next monitoring observation, on MJD 54691.9, 
the total 2.5-25 keV flux had risen to the outburst 
maximum of (6.4 ± 0.8) x 10 -10 erg cm" 2 s _1 , im- 
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plying a peak flux from IGR J00291+5934 of 5.7 x 
10' 10 erg cm -2 s _1 . The daily averages from the Swift 
Burst Alert Telescope (BAT) suggest that this rise oc- 
cured in a single day0 although their large uncertain- 
ties (~3 x 10 -10 erg cm~ 2 s' 1 ) limit the significance 
of this measurement. After «3 d near its maximum, 
the first outburst decayed with an e-folding time of 
(1.8 ± 0.3) d, returning to quiescence in «5 d. The 2.5- 
25 keV fluence of the first outburst, after subtracting 
a 0.68 x 10' 10 erg cm' 2 s' 1 mean contribution from 
V709 Cas, was (3.0 ±0.1) x 10' 4 erg cm- 2 . Most of the 
uncertainty in this figure reflects the uncertainty in the 
V709 Cas flux. 

Following the first outburst, IGR J00291+5934 re- 
turned to quiescence for 30 d. The mean 2.5-25 keV 
flux from the PCA field centered on IGR J00291+5934 
was 0.70 x 10~ 10 erg cm -2 s' 1 during this period, with 
an rms scatter of 0.08 x 10 -10 erg cm' 2 s' 1 . This flux 
and scatter are consistent with the monitoring observa- 
tions before and after the 2008 outbursts and with the 
expected V709 Cas contribution. 

Unlike the first outburst, the second had an approx- 
imately symmetric light curve. In «7 d it rose from 
quiescence to a peak 2.5-25 keV source flux of (4.3 ± 
0.2) x 10' 10 erg cm' 2 s -1 , remained at this peak flux 
for 2-3 d, then dimmed back to quiescent levels over the 
subsequent ~7 d. Observations during the decay were 
too sparse to determine whether the decay was exponen- 
tial. Despite its qualitatively different light curve, the 
second outburst had approximately the same fluence as 
the first: (3.1 ± 0.1) x 10' 4 erg cm" 2 . 

A critical flux of / x , c = 4.4 x 10' 10 erg cm' 2 s' 1 (2.5- 
25 keV) played an important role in the three outbursts 
observed from IGR J00291+5934. The top panels of Fig- 
ure[2]show their light curves, with f xc marked by the dot- 
ted horizontal line. The 2004 outburst was the brightest, 
longest, and most frequently observed by RXTE. While 
its flux was above f x ,c, it decayed linearly at a rate of 
(-0.89 ±0.01) x 10' 10 erg cm' 2 s' 1 d' 1 . Upon reaching 
f x . c , the decay quickened and became exponential, with 
an e-folding time of (1.9 ± 0.2) d (in rea sonable agree- 
ment w ith the 2.2 d time scale reported bv lFalanga et al.l 
I2005bl in their analysis of the INTEGRAL light curve of 
this outburst). The critical flux played a similar role in 
the 2008a outburst: the slow decay seen during the first 
3 d of the outburst quickened considerably after cross- 
ing below f xc . Finally, the symmetric light curve of the 
2008b outburst slowly rose to a maximum flux of ~f x ,c, 
then reversed and slowly decayed. The implications of 
this critical flux will be considered in Section |U 

The linear decay that we report during the first 7 d of 
the 2004 outburst is contrary to theoretical expectations, 
and it therefor e merits close exam ination. The accretion 
disk model of iPowell et al.l (|2007f ). which we discuss in 
detail in Section 14.21 predicts an exponential decay for 
this part of the outburst, and Fala nga et al.l (|2005bl ) re- 
ported an exponential fit to these data. To parameterize 
the shape of this decay, we use the following model: 

u = [fj + ■ e-J ■ (i) 

8 BAT transient monitoring provided by the Swift /BAT team at 
http: / /swif t .gsfc .nasa. gov/docs/swift /results/transients/ 
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Fig. 3. — Light curve for the first 7 d of the 2004 outburst and a 
comparison the best-fit linear (solid line) and exponential (dashed 
line) models. The dotted line at bottom indicates fx,c- 

We fit the initial flux J x q at time t = 0, the decay time 
scale r, and the shape parameter £. This choice of model 
will be motivated by our investation of irradiated accre- 
tion disks in Section [4. 2 [ but for now note that it encom- 
passes linear (£ = —1), exponential (the limit £—*•()), 
and 1/t decay (£ = 1). 

The brightest portion of the 2004 light curve is consis- 
tent with linear decay: £ = -1.05 ± 0.13 with x 2 = 38.2 
for 16 degrees of freedom. In comparison, forcing an 
exponential model (i.e., holding £ fixed at very nearly 
zero) gives x 2 = 121.4 for 17 degrees of freedom. Fig- 
ure [3] compares these two fits. Some irradiated accre- 
tion disk models have £ = 1/4, which fares even worse: 
X 2 = 177.0 for 17 degrees of freedom. The time scale for 
the linear fit is r = 4.9 ± 0.1 d, which is the length of 
time needed for the flux to fall by the amount f x , c - This 
t is similar to the 6.6 d e-folding time scale reported by 
iFalanga et al.l (|2005bf ) for an exponential fit. Note that 
we can exclude an exponential decay even without the 
outlying first observation (on MJD 53342), since the % 2 
for an exponential fit remains roughly three times the 
linear \ 2 without this point. 

3.2. Spectral Evolution during the Outbursts 

The unusual light curves of the 2008 outbursts 
of IGR J00291±5934 naturally raise the question of 
whether its spectrum undergoes significant change dur- 
ing the outbursts. In particular, does the spectrum show 
any transition accompanying the change in the decay rate 

&1 1 fx.c? 

Other sources with "knees" in their light curves do 
not have associated abrupt spectral changes. RXTE 
observations of AMSPs are generally well fit by an 
absorbed ~1 keV blackbody plus power law. In the 
case of SAX J1808.4— 3658, the power law index grad- 
ually softe ned from 1.5 to 1.9 over the courses of the 
outbu rsts ( Gicrlin ski et al.ll2002t llbragimov fc Poutane"nl 
2009), but no significant change in this trend coincided 
with that source's tra nsition from slow to rapid decay 
(|Hartman et al.ll2009bD . For SWIFT J175 6.9-2508, the 
photon index softened from 1.8 to 2.0 (jLinares et al.l 
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Fig. 4. — The spectral fits for the three outbursts. The 2.5- 
25 keV flux from IGR J00291+5934 is shown for reference, in units 
of 10 — 10 erg cm -2 s — 1 . Dotted gray lines mark the critical flux 
f XlC and the times at which the outbursts cross this flux. The 
lower plot shows the fit components: a predominant power law 
and an intermittently detected ~1 keV blackbody and 6.4 keV 
iron line. The blackbody and iron line components are marked 
with circles where they are detected with at least 90% confidence; 
accompanying vertical bars show the times of observations, with 
the width of the bars showing their lengths. 

l2008f) . XTE J0929-314 and XTE J1751-305 had 
roughly constant spectral parameters , with no change 
across the knees of their ligh t curves (jJuett et alj|2003t 
iGierlihski fc Poutanenll2005h . A nalysis of th e 2004 out- 
burst of IGR J00291+5934 by IPaizis et all (p505h re- 
ported an unchanging spectrum: their fits used a ~1 keV 
blackbody component and a constant photon index of 
1.6. They also found hints of a 6.4 keV iron line. 

To measure the spectrum of IGR J00291+5934, we 
first fit the contribution from the nearby intermediate 
polar V709 Cas. For the 2004 outburst, we used the 
RXTE observations of the IGR J00291+5934 field fol- 
lowing its quiescence, and for 2008 we used the obser- 
vations during the 30 d quiescence between outbursts. 
In both cases we fit V709 Cas with an absorbed thermal 
bremsstrahlung spectrum plus a narrow iron line fixed at 
6.4 keV. The resulting plasma temperatures of 26 ±5 keV 
in 2004 and 31 ±4 keV in 2008 agree with earlie r findings 
(|de Martino et al.l 120011: iFalanga etail l2005aD . When 
performing subsequent fits of IGR J00291+5934, we held 
the parameters of the V709 Cas contribution fixed for 
each epoch, effectively subtracting its contribution from 
the total flux. To reduce the uncertainty introduced by 
potential errors in our V709 Cas fits, we only measured 
the spectral parameters of IGR J00291+5934 using ob- 
servations in which it was brighter than V709 Cas. 

We fit the IGR J00291+5934 flux using an absorbed 
blackbody plus power law, with the addition of a 6.4 keV 
iron line when needed. The absorption column was held 
fixed at rig = 4.6 x 10 21 cm -2 based on the optical 
work of iTorres et a l. ( 2008). Figure 2] shows the key re- 
sults. The blackbody temperature was consistent with 
a constant 0.9 keV . In agreement with the analysis of 
IPaizis et al.l (|2005[ ). we found that this blackbody com- 
ponent fades more rapidly than the power law compo- 
nent during 2004, becoming insignificant around the time 
when the total flux falls below f XtC ; however, this is not 
the case for the 2008 outbursts, for which the ratio of the 
blackbody and power law contributions remains roughly 



constant throughout. We found that the power law com- 
ponent softened during the 2004 and 2008a outbursts. 
Data for the 2008b outburst are limited but are consis- 
tent with a constant power law index. This s oftening dur- 
ing 20 04 is contrary to the results reported in lPaizis et al.l 
( 2005) . This discrepancy was likely due to an insufficient 
compensation for the harder spectrum of V709 Cas: if 
its contribution is not subtracted, the overall power law 
index is indeed roughly constant. These 2008 spectra 
are generally compatible with the two contemporane- 
ous Swift XRT detections, which yielded power law in- 
dices of ss 1.6 and did not require blackbody components 
(|Lewis et al.ll2010D . 

No abrupt changes in the spectrum coincide with the 
break in the light curve decay at f X)C . The rate at which 
the photon index gradually increases does not change at 
this point, and no other spectral parameter is signifi- 
cantly affected by this transition. 

3.3. Pulse Timing during the 2008 Outbursts 

The pulse TOAs during the 2008 outbursts exhibit 
a low level of timing noise relative to other AMSPs. 
The rms amplitude of the phase residuals of the best- 
fit constant-frequency model is 0.016 rotational cycles 
(27 fis). Counting noise accounts for approximately 
half of these residuals, contributing an rms amplitude 
of 0.011 cycles. The remaining intrinsic timing noise 
does not exhibit any long-period tr ends or abrupt phase 
shifts (as in SAX J 1808.4— 3658; IBurderi et all [2005 
IHartman et al.l 120081) or correlations wi t h the flux (as 
in XT E J1814-338; iPapitto et all [20071 : iPatruno et all 
2009b). This greatly simplifies analysis, allowing us to 
estimate the uncertainties of timing model parameters 
by scaling the phase uncertainties such that the reduced 
\ 2 of the model is unity. (Because the initial phase un- 
certainties assume only counting noise, which accounts 
for roughly half the total noise, the scaling factors were 
consistently «v2-) A separate bootstrapping analysis 
supports the validity of this approach by producing sim- 
ilar model uncertainties. 

We first developed separate timing models for each of 
the 2008 outbursts. These models and the amplitudes 
of their phase residuals are summarized in Table [TJ with 
a constant-^ model and a non-zero v model fit for each 
outburst. Including a frequency derivative yields a spin- 
up for both outbursts, but at low significance: 1.0 a for 
2008a and 1.8 a for 2008b. The large uncertainties on 
these frequency derivatives are due to the sparsity of ob- 
servations, particularly during the beginnings and ends of 
the outbursts when phase measurements have the great- 
est leverage for fitting ;>, and the decreased effective area 
of the PCA. While the 2008 data are adequately fit by 
a constant-^ model, we cannot conclusively rule out the 
presence of a spin-up of the magnitude o bserved during 
the 20 04 outburst, 8.4(6) x 10~ 13 Hz s _1 (jFalanga et alj 
I2005bn . 

The proximity of the 2008 outbursts enables phase con- 
nection between two AMSP outbursts for the first time. 
With frequency accuracies of ^30 nHz for the two out- 
bursts, the phase uncertainty when connecting across the 
30 d interim of quiescence is ~0.05 cycles — small enough 
to ensure that we do not miss or include an extra rota- 
tion of the neutron star. This long baseline greatly im- 
proves the accuracy of our pulse timing model, and the 
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TABLE 1 TABLE 2 

Frequency models for the 2008 outbursts Orbital parameters for ICR J00291+5934 



v v Residuals 





(Hz) 


(10" 13 Hz s~ 


1 ) (rms cycles) 


2008a 


598.89213081(4) 




0.0148 




598.89213081(4) 


10(10) 


0.0133 


2008b 


598.89213084(3) 




0.0175 




598.89213084(3) 


4.5(2.5) 


0.0150 


Both 


598.892130804(2) 




0.0163 


Note. — For models with 
the middle of the outburst. 


a non-zero v 


v is specified in 



best constant-frequency fit is given in Table [TJ This fre- 
quency is consistent with the frequencies measured for 
the individual outbursts. The addition of a frequency 
derivative does not significantly improve the fit. 

This phase connection between outbursts does require 
the assumption that torques during the 30 d quies- 
cence are small. A quiescent spin derivative of \v\ > 
3 x 10~ 13 Hz s _1 would introduce an unknown number 
of phase wraps. We will address whether such torques 
are physical in the discussion. It is worth noting that the 
phases of the two 2008 outbursts line up quite well. If 
we extrapolate the phase solution of the 2008a constant- 
frequency model forward and apply it to the 2008b out- 
burst without any additional fitting, the residuals during 
2008b have a mean offset of only +0.02 cycles relative to 
our solution. If phase wraps occurred between 2008a and 
2008b, then this excellent alignment would be purely co- 
incidental. 

Finally, we consider the possibility that the marginal 
spin-ups during the 2008 outbursts are real. For con- 
cision, we compare their effects relative to the best-fit 
constant frequency during 2008, v = 598.892130804 Hz. 
The values of v listed in Table [1] would result in a fre- 
quency of v = v$ + 0.23(13) /iHz at the end of 2008a 
(at MJD 54697.0) and v = v - 0.24(11) ^Hz at the 
beginning of 2008b (at MJD 54730.0). The neutron 
star would have to undergo a mean spin-down of i> = 
— 1.7(6) x 10 -13 Hz s _1 during the intervening 33 days 
of quiescence to account for this decrease of frequency. If 
a spin-up is present during these outbursts, then a spin- 
down of a similar magnitude necessarily follows them. 

3.4. Long-term Spin and Orbital Evolution 

Despite the highly accurate (2 nHz uncertainty) tim- 
ing solution attained by phase connecting the 2008 out- 
bursts, it is still not feasible to phase connect the 2004 
and 2008a outbursts. During the 3.6 years of quiescence 
separating these outbursts, magnetic dipole torques from 
the neutron star's unknown field strength (say, 10 7 ' 5 - 
10 8 - 5 G) will produce 1-50 phase wraps relative to a 
model with a constant frequency. This prevents the con- 
struction of a coherent timing model for all the data. 

Comparing the frequencies of the outbursts is mor e 
useful. From the ephemeris of iFalanga et al.l j2Q05b), 
the frequency at the end of the 2 004 outburst will be 
v = v Q + 0.32(4) jitHz; the model of iBurderi et~all (|2007l ) 
and our own analysis give similar figures. As discussed 
in the previous section, the determination of v for the 
2008a outburst is more ambiguous. If we accept the 



Orbital period, P olb (s) 8844.076729(9) 
Orbital period derivative, P orb (95% UL; s s _1 ) < 3 X 10" 11 

Projected semimajor axis, a x sini (light- ms) 64.9911(4) 
Time of ascending node, T asc (MJD, TDB) 54000.0739971(7) 
Eccentricity, e (95% confidence upper limit) < 2 X 10~ 4 

marginally significant spin-up, then the frequency at the 
beginning of the 2008a outburst is v — vq — 0.21(19) /iHz; 
if not, then we can phase connect the 2008 outbursts to 
give the constant frequency of v§ with an uncertainty 
of 0.002 fjRz. These scenarios produce respective fre- 
quency drops of 0.53(20) yuHz and 0.32(4) /iHz during the 
3.6 years of quiescence. Averaged over the entire quies- 
cent period, these drops represent a mean spin-down of 
-(3-5) x 10~ 15 Hz s" 1 . 

Connecting the orbital phases of the 2004 and 2008 
outbursts provides a great increase in the accuracy of the 
orbital period. Table [2] summarizes the orbital ephemeris 
derived by combining the data from all three outbursts. 
Detecting an orbital period derivative is possible in prin- 
ciple with three orbital phase measurements, but the 
proximity of the 2008 outbursts prevents anything more 
than a weak upper limit on the orbital period derivative. 
The effect of any physically plausible derivative on the 
orbital phase during the 30 d separating the outbursts 
would be too small to measure. 

3.5. Pulse Profiles 

The pulses from IGR J00291+5934 were entirely con- 
sistent with a sinusoid across the 2-60 keV energy band 
of the PCA. A relatively low amount of timing noise al- 
lowed us to integrate the pulse profile over both 2008 
outbursts. Folding all the 2.5-15 keV photons gave a 
profile that was well-fit by a pure sinusoid with no har- 
monics (x 2 = 259.8, 253 degrees of freedom). The frac- 
tional amplitude in this band was 10% rms and remained 
roughly constant throughout both outbursts. The 95% 
upper limits on the fractional amplitude for the second 
and third harmonics were both 0.36% rms, or <4% of 
the fundamental's amplitude. These results are consis- 
tent with the non-detection of harmonics during the 2 004 
outburst (jGallowav et al.ll2005t IFalanga et al.ll2005bl ). 

To measure the energy dependence of the pulsations, 
we divided the PCA response into 15 energy bands and 
folded all the 2008 outburst data for each. Figure[S]shows 
the measured pulse arrival times and fractional ampli- 
tudes. Data from the 2004 outburst are also shown for 
comparison. Pulsations were detected with 99% confi- 
dence at 2-35 keV (2-48 keV during the 2004 outburst). 

From 2-8 keV, the pulse arrival times follow the usual 
pattern of soft lags that is se en in most other AMSPs. As 
IFalanga k. Titarchukl (|2007f ) noted for the 2004 outburst, 
above 8 keV the trend reverses and harder bands increas- 
ingly lag. Agreement between the 2004 and 2008 phases 
is generally good. There was no significant change in the 
magnitude of the lags over the course of the outbursts, 
in contrast to the flux dependence of the lag s observed 
in SAX J1808.4-3658 (|Hartman et al.ll2009bji . 

The fractional amplitude decreases significantly until 
8 keV, above which the energy dependence is weaker but 
with some evidence of further decrease. The error bars 
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5 10 20 

Energy (keV) 

Fig. 5. — The phases and fractional amplitudes as a function of 
energy. The black points are integrated over both 2008 outbursts; 
the gray points are for the 2004 outburst. Error bars show the 
statistical uncertainties. For the fractional amplitudes, an addi- 
tional uncertainty in the overall vertical offset is introduced by the 
approximated flux from V709 Cas. 

for both the phases and amplitudes only account for sta- 
tistical uncertainty. The fractional amplitude errors do 
not include the uncertainty in the flux from V709 Cas, 
which is assumed to be the same for both outbursts. 
Nevertheless, the offset between the 2004 and 2008 am- 
plitudes is significant. The V709 Cas flux would have 
to be 30% higher during 2004 to account for this frac- 
tional amplitude difference, but a comparison of its flux 
from monitoring observations after the 2004 and 2008 
outbursts returned to quiescence show that this was not 
the case. 

No pulsations were detected during the quiescent pe- 
riod between the 2008 outbursts. Folding all data in 
MJD 54700-54730 yields a 95% upper limit on the pulse 
fractional amplitude of 0.4% rms, accounting for instru- 
mental background only. This limit further supports the 
assumption that the non-background photons during this 
period are from V709 Cas, and IGR J00291+5934 is in- 
deed in quiescence. 

3.6. Aperiodic Timing 

A broadband power spectrum of IGR J00291+5934 re- 
veals a number of notable features. Figure |5] shows the 
spectra, and Table [3] lists some of their properties. For 
comparison, the top plot of Figure [6] shows the power 
spectrum of the first 6 d of the 2004 out burst, a selec- 
tion la beled "subset Al" in the analysis of iLinares et al.l 
(2007). High levels of aperiodic timing noise distinguish 
the spectra of this source from those of other NS LMXBs. 
The noise level is flat down to very low frequencies 
(^0.01 Hz), resulting in integrated fractional variabili- 
ties over 0.01-100 Hz of 40-60% rms. Additionally, two 
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Fig. 6. — The timing spectra of the three outbursts from 
IGR J00291+5934. Dotted gray curves show Lorentzian compo- 
nents centered at the origin (i.e., with a coherence of Q = 0); solid 
gray curves show QPO components, modeled with Lorentzians with 
non-zero coherences. Subset Al of the 2004 outburst includes only 
the b rightest part of that outburst, as defined by Linar es" et al.l 
(2007)/ 



harmonically related quasi-periodic oscilla tions (QPOs) 
are p resent at around 0.02 Hz and 0.04 Hz (jLinares et al.l 
1200% . 

The broadband shape and integrated variability of the 
2008 power spectra are similar. The spectra are roughly 
flat down to a break at ~0.01 Hz, and the integrated 
fractional variability is ~50% rms. No QPOs were de- 
tected during the 2008a outburst, but its higher overall 
noise level would be sufficient to bury the QPOs observed 
during 2004 if they were present at the same fractional 
amplitude and coherency («5% rms; Q ~ 5), so their 
non-detection is not constraining. During the 2008b out- 
burst, a single QPO was definitively seen near the low- 
frequency break. It is unclear which harmonic it repre- 
sents if it is a member of a harmonically related QPO 
pair: at 0.026 Hz, it falls between the QPO frequen- 
cies seen in the earlier outburst. An additional QPO at 
0.45 Hz is needed to model the steeper power spectrum 
"hump" present during 2008b. This QPO has an ampli- 
tude of (13 ± 2)%, a coherency of Q = 1.5 ± 0.5, and 
is present in all 2008b observations with sufficient inte- 
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TABLE 3 

Aperiodic timing properties of IGR J00291+5934 



2004 ; 



2008a 



2008b 



Interval of spectrum (MJD) 


53341.0-53348.2 


54694.0-54695.0 


54733.0-54737.5 


Min/mean/max flux (10 -10 erg cm -2 s — 1 ) 


5.4 / 6.8 / 10.9 


4.0 / 4.2 / 4.8 


3.5 / 4.2 / 4.8 


Total 0.01-100 Hz variability (% rms) 


38.8 ±0.2 


50.5 ±0.5 


51.3 ± 0.7 


Low-frequency break (10 — 2 Hz) b 


4.9 ±0.2 


2.1 ±0.2 


3.3 ± 0.2 


Lower QPO (10~ 2 Hz) 


2.03 ±0.06 




2.64 ±0.05 


Upper QPO (10~ 2 Hz) c 


4.37 ± 0.08 




45 ± 2 



a Subset Al only. The 2004 figures are from our fit but agree with Linares ct al. (2007). 

b The break frequency given is !^ max (sec Bclloni ct al. 2002) of the lowest- frequency zero-centered 

Lorentzian. 

c Unlike the 2004 QPOs, the QPOs during 2008b arc not related. 



gration time to detect it. The 20 08 data suppo r t som e 
of the overall trends identified by ILinares et al.l (|2007l ). 
particularly the anticorrelation between the overall vari- 
ability and the flux and the positive correlation between 
the break frequency and the flux. 

4. MODELING THE OUTBURST LIGHT CURVES 

The most unusual feature of the 2008 activity of 
IGR J00291+5934 is its light curve: the source under- 
goes two outbursts in rapid succession. We are unaware 
of any other NS LMXB for which this is the case0 This 
double outburst provides a rigorous test of the accretion 
disk models put forward to explain the transient nature 
of many LMXBs. A successful theory must accomplish 
the following: 

Explain the light curves of the 2004 an d 2008a out- 
bursts, particularly the knee at the critical flux f x , c - 
These two outbursts were preceded by a long (~3 yr) 
period of quiescence, and they followed a fast rise / slow 
decay profile. Both showed a knee in the decay at a crit- 
ical 2.5-25 keV flux of / X)C = 4.4 x 10~ 10 erg cm" 2 s _1 . 
During the 2004 outburst, the decay prior to the knee 
was nearly linear, with a timescale of f x ,c/ fx — 4.9 d; 
during the 2008a outburst, a lower peak flux and sparser 
sampling prevented characterization of the nature of the 
pre-knee decay. Both outbursts showed exponential de- 
cays after the knee with e- folding times of «1.9 d. 

Stop the 2008a outburst before the disk becomes too de- 
pleted to fuel the 2008b outburst. Unless the mass trans- 
fer rate from the companion can vary by orders of mag- 
nitude on a ~1 month time scale, the mass needed for 
2008b could not have accumulated during the 30 d qui- 
escence that preceded it. The 2.5-25 keV fluence of the 
2004 outburst was 5.2 x 10~ 4 erg cm~ 2 s _1 , and the 
summed fluence of the 2008 outbursts was similar, at 
6.1 x 10 -4 erg cm~ 2 s _1 . Taking these fluences and a 
3 yr recurrence time as typical, we estimate that at least 
95% of the matter accreted during 2008b must have been 
present in the disk at the end of 2008a. 

Account for the symmetric light curve of 2008b, and 
why its peak flux is f x ,c- In contrast to outbursts fol- 

9 Some NS LMXBs do show much shorter (< 3 d) and dim- 
mer mini-outbursts or fl ares that recur frequentl y; XTE J1751— 305 
UMarkwardt et alj|200?tl and NGC 6440 X-2 JHeinke et al.ll20ToT) 
are notable examples. At the other extreme, the intermittently 
pulsating AMSP HETE J1900.1-2455 has been in outburst since 
2005, with the e xception of 1—6 d in 2007 during which it was 
briefly quiescent (Dcgenaar ct al. 2007). Unlike these systems, all 
the outbursts from IGR J00291+5934 resemble the typical ~10 d 
outbursts seen in most NS LMXB transients. 



lowing a long period of quiescence, the 2008b outburst 
rose slowly over «7 d until reaching a maximum flux 
of f XtC . After 2-3 d at this maximum flux, the outburst 
then slowly decayed back to quiescent levels over another 
«7d. ' ' 

We consider this unique phenomenolog y in light of 
three models: the disk diffusion model of IWood et al.l 
(2001), which was developed to explain a similar double 
outburst from the black hole transient XTE J1118+480; 
accretion from a partially ionized irradiated disk; and ac- 
cretion at rates near the onset of a quasi-propeller state, 
in which the centrifugal acceleration of infalling matter 
by the NS magnetosphere inhibits but does not halt ac- 
cretion. 

4.1. Disk Diffusion Model 

The one other LMXB with a published report of a dou- 
ble outburst is the black hole transient XTE J1118+480. 
Its first outburst lasted 40 d and followed a fast rise / ex- 
ponential decay profile. After 30 d quiescence, a second 
outburst lasted 150 d and had an irregular profile with 
multiple peaks (jWood et al.ll200lD . 

Additionally, IGR J00291+5934 and X TE Jl 118+480 
have s imilar aperiodic timing properties. ILinares et al.l 
(|2007t) note that the strong very-low-frequency (< 
0.1 Hz) variability and high overall variability of 
IGR J00291+5934 much more closely resemble the tim- 
ing spectra of black hole LMXBs than other NS LMXBs. 
In section 13.61 we reported a 0.45 Hz QPO during 
the 2008b outburst of IGR J00291+5934; scaled for the 
difference in central object masses (1.4 M Q vs. 6- 
7 Mq), it matches the ~0.1 Hz QPO obs erved during 
the s econd outburst of XTE J1118+480 (|Wood et all 
120001) . The data for the 2008b outburst were not ad- 
equate to determine whether the QPO frequency in- 
creased over the course of the outburst, as was observed 
for XTE J1118+480. 

Given the observational similarities between these 
so urces, it is w o rth co nsidering the disk diffusion model 
of IWood et~aLl (|200lD . Noting that the disk instability 
model was insufficient to explain the double outburst of 
XTE Jl 118+480, they proposed that the light curve fol- 
lowed from a varying mass loss rate from the companion, 
which would accrete onto the NS after diffusing through 
the accretion disk. The resulting light curve would be 
the convolution of the companion's mass loss function 
and a fast-rise, exponential-decay diffusion function that 
acts on the viscous time scale of the disk. It is un- 
likely that this model can account for the outbursts of 
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IGR J00291+5934, however. If the viscous time scale of 
the disk is similar to the decay time scales of the 2004 and 
2008a outbursts, then their observed light curves would 
require a very rapid (<3 d) transfer from the companion 
to the outer accretion disk of most of the mass consumed 
in those outbursts. Furthermore, a disk diffusion model 
cannot explain the knee in the light curves. The 2008b 
outburst would require only slightly slower transfer of 
matter from the companion to the disk. 

4.2. Accretion from an Irradiated Disk 

We next consider the outbursts from the perspective 
of the disk i nstability mod e l as a pplied to soft X-ray 
transients by [King k Ritterl (jl998l hereafter KR98) and 
iPowell et all (|2007L hereafter P07). In this approach, the 
accretion disk is divided into two regions: a hot inner 
disk in which ionization results in a high viscosity and 
the inward migration of matter; and a cooler, non-ionized 
outer disk in which the inward drift is much lower. 

The radial extent of the hot, ionized region can be 
calculated by considering the temperature of the disk 
due to irradiation from the NS. If an accretion disk with 
a scale height H oc R n and albedo 77* is i rradiated by a 
small central source with a luminosity L^. lde Jong et aLl 
(|1996f) calculate its temperature profile to be 



rpA 



1 



H 



AttR 2 R 



(n - l)L x . 



(2) 



Following the argument of P07, the hot region will extend 
out to some temperature at which the gas in the disk 
becomes mostly ionized. If we denote the radius of this 
temperature as Rh, we can solve to find 



R 3 h ~ n ex L x 



oc M x . 



(3) 



where M x is the accretion rate onto the compact ob- 
ject. When the disk is illuminated by a point-like central 
source, as is the case here, n « 9/7 (KR98). 

Above some critical accretion rate M XyC , the luminosity 
will be sufficient to ionize the entire disk. Thus we can 
write 

/ . x l/(3-n) 

"" (4) 



Rd 



isk 



To simplify, we define the dimensionless accretion rate 
rh x = M x /M x ^ c and the disk opening-angle parameter 

7 = 1/(3 — n). As M x falls below M X]C , the transition 
from a fully to a partially ionized disk will generally cause 
the decay rate to increase, causing a knee in the light 
curve. This mechanism has been invoked to explain the 
knee seen in the outburst light curves of many transient 
LMXBs (e.g., P07). 

4.2.1. Assumptions of the Viscosity Model 

The models of KR98 and P07 make two assumptions 
about the disk viscosity. (1) The viscous time scale of 
the ionized region of the accretion disk is shorter than 
the decay time scale of the outburst: 



Rl 



< 



d\nM x 



dt 



(5) 



Here v is some mean kinematic viscosity of the ionized 
region. (Throughout this discussion we adopt the con- 
vention that v represents the kinematic viscosity rather 
than the NS spin frequency, unless otherwise noted.) If 
this condition is met, the surface density of the hot disk 
will relax into the quasi-steady profile 



Jhot 



(R) 



M x 
Sttv 



(6) 



(2) The viscous time scale of the cold region is far longer 
than the decay time scale, so its surface density profile 
can be treated as independent of M x . 

For simplicity, the calculations of these papers also 
make the approximation n = 1 (or equivalently 7 = 1/2) 
for a NS LMXB, and they use a constant mean kine- 
matic viscosity for v. We adopt a more general approach 
by keeping 7 a free parameter and by allowing the disk 
viscosity to vary with the radius and accretion rate: 

v(R,M x ) = v {R/R A ^ymi. (7) 

For the IShakura k Svunvaevi (|1973l ) disk solution, £ = 
3/10 and /? = —3/4. For an irradiation-dominated disk, 
we assume v oc T, and eq. ([2]) gives £ = 1/4 and (3 = 
(n-3)/4« -3/7. 

4.2.2. Light Curve Decay of a Fully Ionized Disk 

Following the arguments of KR98, we can analytically 
solve for the light curve when the accretion disk is en- 
tirely ionized: 



M x = M 



m: 



t 

' Thot 



-i/S 



(8) 



The initial accretion rate is M x ,o a ^ time t = 0, and the 
time scale for decay is 



Thot 



2(1 



R 2 



disk 



3(2-/3) is 



(9) 



Note that in the limit of a uniform and constant viscosity 
(/3 = £ = 0), these equations reduce to the exponential 
decay given in KR98: M x = M xfi exp(-t ■ 3z//i^ isk ). 

In this model, a knee in the decay of the outburst is 
generally attributed to a transition from a fully to par- 
tially ionized disk. However, it is not possible to reconcile 
this hypothesis with the linear decay observed when f x > 
f X)C during the 2004 outburst of IGR J00291+5934. Fit- 
ting eq. (j8]) to the light curve, we find £ = — 1.05±0.13, in 
strong disagreement with the expected value of £ = 0.25 
for an irradiated disk. Therefore the accretion disk of 
IGR J00291+5934 is not fully ionized during the initial 
decay of the 2004 and 2008a outbursts, and the signifi- 
cance of f XiC must be otherwise explained. 

4.2.3. Behavior of a Partially Ionized Disk 

When the disk is only partially ionized, we consider 
the mass budget of the hot disk, following P07: 



Mhot = —M x + lJ, C oid(Rh) + 2wRhRhZ(Rh) 



(10) 



The Mhot term can be derived by integrating the hot-disk 
surface density out to Rh, then taking the time deriva- 
tive. From assumption 1, this surface density follows 
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eq. ©. The Hcoid(Rh) term gives the viscously driven 
inward flow of matter from the cold disk at Rh ■ From as- 
sumption 2, the radial flow in the cold disk is far smaller 
than in the hot disk, so we can neglect this term. The 
'E(Rh) term reflects the addition of cold-disk matter due 
to the encroaching ionization radius during the outburst 
rise, and conversely the return of matter to the cold disk 
as Rh recedes during the outburst decay. 

Applying these relations to eq. (fTO)) and casting it in 
terms of the dimensionless accretion rate m x = M x /M XtC 
gives a general model for the partially ionized disk: 



27T7m 



h 2 T 



M. r r 



1 



C • -a- 



-i -l 



(11) 

To simplify, we define the constant ( = £ — j(2 — /?). 
The behavior is determined by the form of H(Rh). For 
Rh < 0, the surface density of the hot disk is relevant: 
E(i?) = Ehot(-R) from eq. ©. For R h > 0, E(i?) is the 
surface density of the cold disk, which depends on the 
state of the disk prior to outburst. As the forms of Ehot 
and E co id are generally different, we must handle the two 
cases separately. 

During the decay of an outburst from a partially ion- 
ized disk, the fall of M x will cause Rh to move inward. 
As it does, the surface density of matter that this cooling 
front encounters will follow the quasi-stable hot disk pro- 
file of eq. (|6]), and from assumption 2 this ^hot(Rh) wm 
become "frozen" into the cold disk profile as Rh passes 
inward. 

This "freezing in" of the hot disk surface density pro- 
vides another test of the hypothesis that the light curve 
knee marks the transition from a fully to partially ionized 
disk. The resulting density profile left behind is 

-/3+U-0/7 



E(i*) = 



R 



\i?disk 

For the 7, £, and ft of an irradiation-dominated disk, 

ThotMx.c ( R 



(12) 



E(i?) = 0.5 x 



fi> 2 

-"■disk 



(13) 



, -Rdisk , 

If we assume full ionization at / x ,o then we can measure 
Thot = 6.6 d from the initial decay of the 2004 outburst 
(jFalanga et all l2005bf ) and use f x /f x , c = M x /M XiC to 
convert between flux and the relative accretion rate. In- 
tegrating eq. (|13l) over the disk then gives a post-outburst 
mass sufficient to account for only 70% of the observed 
fluence during the 2008b outburst. No physical viscosity 
model will leave 95% of the 2008b fluence remaining in 
the disk after the 2008a outburst, as required. Again, it 
is difficult to reconcile the observed outbursts with the 
accretion disk of IGR J00291+5934 ever being fully ion- 
ized. 

In contrast, this approach can successfully model the 
outburst light curve when f x > f x>c if we assume a par- 
tially ionized disk. In this scenario eq. (fTTj) can be solved 
analytically, giving a decay for a partially ionized disk 
with the same form as the fully ionized case but with £ 
as its shape parameter: 

-i/C 



M x 



M 



< 



M~i ■ C— 

Ttiot 



(14) 



Assuming a constant viscosity and a disk scale height 
that increases linearly with radius, we obtain £ = — 1, 
giving the linear decay predicted by KR98 and P07. 
For an irr adiation-dominated di s k and a gas-pressure- 
dominated [Shajmra&IS^un^ae3 (|1973l ) disk, the shape 
parameter is C = —1-17 and ( = —1.30, respectively. The 
observed decay shape parameter of —1.05 ± 0.13 when 
fx > fx,c is compatible with the expected irradiated disk. 
However, another mechanism must be invoked to explain 
the knee at f x . c . 

4.3. Magneto spheric Inhibition of Accretion 

The drop in the accretion rate when the flux falls 
below f XjC may also be caused by the onset of mag- 
netospheric effects. This explanation has been in- 
voked to explain similar drops i n the light curves of 
other NS LMXBs (e .g., Aql X-l. iCampana et alJU^9l 
SAX J1808.4-3658, IGilfanov et all 11998ft . A~sthe ac- 
cretion rate declines, the magnetospheric radius moves 

outward: R m oc M x ' . When R m exceeds the co- 
rotation radius R co , at which the Keplerian orbital fre- 
quency equals the NS spin frequency, infalling matter 
must accelerate to co-rotate with the magnetic field. 
At this point the source is typically said to enter the 
"propeller" regime, as it was originally though that cen- 
trifug al acceleration would eject matter from the sys- 
tem |Illarionov^Su^ae3 [l975|). In fact, this picture 
is not energetically self-consiste nt, and R m /R co > 1.2 
is required for m ass ejection (Rappap ort et al.l 12004 
iPerna et al.l l2006h . Below this limit, matter will build 
up in the vicinity of R m until its material pressure 
pushes the magnetosphere inward, allowing matter to 
accrete onto the NS a nd relieving the pressure on R m 
l|Spruit fc Taaml Il993l ). We refer to this intermediate 
state as the "quasi-propeller" accretion mode, because 
the centrifugal acceleration of infalling matter throttles 
but does not entirely stop accretion. 

Modeling the throttling of the accretion rate induced 
by the onset of a quasi-propeller state is beyond the 
scope of this paper, but we consider some qualitative 
predictions. In particular, it can account for the differ- 
ent shapes of the outburst light curves during 2004 and 
2008a, which both exhibit a fast rise and a linear-then- 
exponential decay, and 2008b, which has a roughly linear 
rise, a plateau, then a roughly linear fall. The key to this 
difference is the distribution of matter in the accretion 
disk prior to the outbursts. Additionally, the rapid drop 
in flux caused by the onset of the quasi-propeller state 
allows enough mass to remain in the disk following the 
2008a outburst to fuel the 2008b outburst. 

First, consider the 2004 and 2008a outbursts. An ac- 
cretion disk that has been quiescent for a sufficient length 
of time will relax into a constant-density state (KR98): 



E(i?) (xH ocR n « R 



1.2 



(15) 



Distributing an initial disk mass accordingly and inte- 
grating eq. flTTl) . we get a light curve that rises faster 
than exponentially. This rise breaks our first assump- 
tion, that the hot disk can maintain a quasi-steady-state 
density profile, but it is qualitatively instructive: it pre- 
dicts the fast rise observed for an outburst following a 
long quiescence. For these outbursts the rising flux never 
became bright enough to fully ionized the accretion disk, 
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causing the initial outburst decay to be approximately 
linear, as discussed in the previous section. Finally, when 
the flux reached f XlC magnetospheric throttling becomes 
effective, causing the decay rate to greatly increase and 
ultimately quenching the outburst. 

The accretion disk that was present at the beginning 
of outburst 2008b would not have had time to relax into 
a quiescent profile, so it depended on the profile left be- 
hind by 2008a. Consider the instantaneous state of the 
accretion disk when the 2008a outburst reached the knee 
a t fx — fx.c- In the cold region of the disk (i.e., outside 
the Rh corresponding to the flux f XjC ), the disk profile 
left behind by the cooling front will follow eq. ([TB")) : 

E(JZ) oc R-^-O/i ~ for R > Rh (f x c ) . (16) 

In the ionized region, the quasi-steady surface density 
profile of cq. (|6|) will be present: 

« M*. oc ir"« R 0A for R<R h (f x , c ). (17) 

As the flux continued to fall, the onset of magnetospheric 
throttling caused the decay time scale to becomes faster 
than the viscous time scale. The result is a greater 
amount of mass remaining in the inner disk than pre- 
dicted by eq. (H3J): 40% more matter would have been 
accreted if the light curve had continued to follow a linear 
decay rather than becoming exponential. The distribu- 
tion of the remaining inner disk mass will fall between its 
state at the onset of magentic throttling and the distri- 
bution left behind by freezing in the quasi-steady profile: 
0.4 < dlogS/dlogi? < 1.7 for R < R h (f x , c ). Finally, 
during the 30 d quiescence the disk profile would begin 
to relax toward the constant-density R distribution, 
but it likely would not have enough time to reach it. 

This disk profile has two repercussions. First, the 
faster shutoff of accretion should leave sufficient mass in 
the disk to fuel the 2008b outburst. Second, the distri- 
bution of this mass can result in the observed 2008b light 
curve. of the inner disk has a smaller gradient than 

it would after a long period of quiescence, so the amount 
of mass that Rh encounters as it expands outward in- 
creases more gradually. As a result, we do not get the 
sudden, super-exponential brightening seen for 2004 and 
2008a. Once the flux reaches / x>c , however, the heating 
front will encounter a lower surface density, as the disk 
ionized by fluxes higher than f XjC will have been more 
fully depleted by 2008a. This change in d log £/d log i? 
will stop the advance of Rh, initiating the outburst de- 
cay. The result is a slow rise, slow decay outburst that 
peaks at f x>c , as observed. 

We have not attempted to explain what stops the rapid 
rise of the 2004 and 2008a outbursts. One possibility is 
self-shadowing by the accretion disk: a warped disk or a 
point at which the radial dependence of the disk height 
decreased would halt or slow the outward movement of 
the ionization front during the outburst rise. Regard- 
less of the mechanism, the lower peak flux of the 2008a 
outburst caused it to more quickly reach the onset of 
rapid decay at f x ,c, ultimately leading to a shorter out- 
burst with an accreted mass of roughly half the mass 
consumed during the 2004 outburst. It is likely that the 
lower peak flux of 2008a was a necessary condition for 
enough matter to be left in the disk to fuel the 2008b 



outburst. 

5. EVIDENCE FOR CHANGES IN THE ACCRETION STATE 

The previous section showed how the transition of 
IGR J00291+5934 from steady accretion to a quasi- 
propeller state at the light curve knee could success- 
fully account for the other features of the outburst light 
curves. This explanation has two requirements of the 
accretion rate corresponding to the knee: first, it must 
not fully ionize the accretion disk; and second, it must be 
compatible with the expected accretion rate at which the 
NS magnetosphere begins transferring angular momen- 
tum to the infalling material. In this section we address 
these requirements. We also consider other evidence sug- 
gesting that IGR J00291+5934 enters a quasi-propeller 
state during its outbursts. Throughout, we compare with 
the 401 Hz AMSP SAX J1808.4-3658, which also shows 
light curve knees that are most likely ass ociated with the 
onset of the quasi-propeller state (e.g., lHartman et al.l 
I2009U rPatruno et al.|[2009ah . 

5.1. The Accretion Rate at the Light Curve Knee 

Di stance estimates to IGR J00291+5934 cover 2-5 kpc 
(see iTorres et al.l 120081 for a recent review); we choose 
4 kpc for our calculations, which is consistent with mass 
transfer arg uments and the lack o f observed thermonu- 
clear bursts ()Gallowav et alJfeOOSf) . The bolometric cor- 
rection factor for our 2.5-25 keV fluxes is Cboi = 2.54 
([Galloway et al.lfeOOSh . With an assumption of a canon- 
ical NS mass and radius (M x = 1.4 Mq, R x = 10 km), 
the accretion rate at the light curve knee is 

4,c = 4 ^y* c hol f x , c = 1.8 x 10- 10 Mq yr- 1 . (18) 

The distance uncertainty introduces a factor of ~2 un- 
certainty in this figure. 

The outburst light curves of SAX J1808.4-3658 also 
show a knee at which decay steepens. Assum ing a dis- 
tance of 3.5 kpc (|Gallowav fc Cumm ing 2006) and iden- 
tical NS parameters, the accreti on rate at that source 's 
knee is also 1.8 x 10~ 10 Mq yr" 1 (lHartman et al.ll2009bl) . 
The perfect agreement of these numbers is a coincidence, 
of course, but it does strongly suggest that the light 
curve knees of the two sources arise from the same phys- 
ical process. Considerable evidence points to this knee 
marking the beginning of the quasi-propeller state for 
SAX J1808.4-3658. A -1 Hz QPO with a very high am- 
plitude (sometimes > 100% rms) is present after the knee 
of some outbursts. iPatruno et al.l (|2009a|) argue that a 
natura l explanation is provided by the ISpruit fc Taaml 
( 1993) instability, wherein matter that accumulates near 
the magnetosphe ric radius quasi - periodi cally drips onto 
the NS surface. lHartman et al.l (|2009bl) show that the 
behavior of the soft lags of the pulsations also changes 
at this critical flux in a way that suggests a change in 
the accretion column geometry that would occur when 

Rm ~ ^co ■ 

Similar evidence suggest that a quasi-propeller state 
turns on in IGR J00291+5934 at fluxes below f XiC . A 
~0.5 Hz QPO was present throughout the 2008b out- 
burst, for which the flux was < f x<c at all times. Al- 
though its 13% rms amplitude was much less than the 
amplitude of the 1 Hz QPO of SAX J1808.4-3658, it 
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is plausible that it too is due to the Spruit-Taam insta- 
bility. The energy dependence of the IGR J00291+5934 
pulsations remained constant through the outbursts, so it 
offered no evidence either way. We will consider other ef- 
fects indicating a strong disk-magnetosphere interaction 
later in this section. 

The hypothesis that the accretion disk of 
IGR J00291+5934 is never fully ionized during the 
observed outbursts is perhaps more controversial. In 
section 14.2.21 we argued that the observed linear decay 
during the first 7 d of the 2004 ruled out a fully ionized 
disk. The same is true for SAX J1808.4-3658: we find 
that the high-flux, slow-decay stages of its 1998, 2005, 
and 2008 outburst light curves are better fit with linear 
models than exponential models^ while the rapid decay 
following the kne es of its outbursts a re approximately 
exponential (e.g., [Gilfanov et al. 1998). Aql X-l shows 
a similar transi t ion t o exponential decay as it fades 
(|Campana et al.l I1998D . On the other hand, the cal- 
culations of KR98 predict that short-period transients 
should fully ionize their disks, and full ionization has 
been invoked to account for the lig ht curves of many NS 
LMXBs (e.g. , IShahbaz et allll99a P07) . 

The 2004 outburst had a peak flux of 2.5 / XjC , implying 
an accretion rate of 4.5 x 10 -10 M Q yr _1 . KR98 predicts 
that the entire disk will be ionized above a critical rate 
of 



a peak-to-peak separation of Av p 



650 ± 40 km s" 



M XiK R = 4.1 x 10 



-10 



Rdi, 



10 6 km 



M e yr" 



(19) 



The disk radius of IGR J00291+5934 must be smaller 
than 10 6 km for the observed orbital parameters and re- 
alistic masses, causing an apparent contradiction. That 
said, the disk model used by KR98 to derive the above 
critical rate is quite basic and likely requires substantial 
modification to produce a physical disk structure in a 
tight binary. A significant shortcoming is the assumption 
of a constant n = d log H/ d log R to describe the vertical 
profile. Simulations of irradiated di sks produce profile s 
that flatten out at larger radii fe.g-. IDubus et al.l ll999). 
screening the outer disk from irradiation and increas- 
ing the M x required to ionized the non-screened portion 
of the disk. Additional ly, irradiation will tend to warp 
LMXB accretion disks (|Pringldll996D . allowing irradia- 
tion of some outer regions while shadowing others. The 
profile of an irradiated disk remains an open question, 
but it most likely can accommodate a partially ionized 
disk even at high accretion rates. 

Finally, we must consider the detection of Ha emis- 
sion during the 2004 an d 2008 outbursts (jTorres et al.l 
120081 : ILewis et"aT] |2010[ ) under the assumption of par- 
tial ionization. The Ha lines were double-peaked and 
showed no variation with orbital phase, clearl y indicat- 
ing their origin in the disk. ITorres et all (|2008l ) reported 

10 The 2002 outburst of SAX J1808.4-3658, which had a peak 
flux 50% brighter than the other outbursts, initially followed an ap- 
proximately exponential decay until reaching a first knee while still 
at a relatively high flux. The decay then slowed, but whether its 
subsequent form was better fit with an exponential or linear model 
depends on the region fit. After a second knee at that source's 
cr itical flux, it entering rapid decay (see Fig. 3 of P07 and Fig. 1 
of IHart man et al. 2008). The peak flux of this outburst possibly 
does fully ionize the accretion disk, with the first knee marking the 
transition to partial ionization. 



on 2004 Dec 5 (MJD 53344), near the peak of the 2004 
outburst. Assuming that this Av P p measured the line- 
of-sight velocity at the edge of a fully ionized accretion 
disk, they estimated an orbital separation of a — (7.2- 
8.3) x 10 5 km, an inclination of i = 22°-32°, and a com- 
panion mass of M 2 = 0.07-0.11 M for a 1.4 M Q NS 
(0.09-0.13 M Q for a 2.0 M NS). If the assumption of 
full ionization is dropped, then their calculated a and 
Mi ranges become lower limits, and their i becomes an 
upper limit. A weaker upper limit on M% is provided 
by the requirement that the com panion fits within its 
Roche lobe, Riving M 2 < 0.25 M (jGallowav et"aL|[2005l: 
ITorres e~a l. 2008). The resulting range of allowed incli- 
nations is 10°-32°. A low inclination for this system is 
supported by th e lack of harmoni c content in the X-ray 
pulsations fe.g. JMuno et al.ll2002fl . 

5.2. The Magneto spheric Radius at the Knee 

As discussed at the beginning of Section 14.31 the mass 
accretion (or ejection) regime depends on the magneto- 
spheric radius, 



R,, 



k m (2GM x )-^ 7 M- 2 / 7 ^ 7 . 



(20) 



Here \i is the magnetic dipole moment, and k m is 
an order-unity constant that encapsulates the complex 
disk/magneto sphere interaction . Magnetohydrodynamic 
simulations of Lon g et al.l (|2005f ) suggest a value of k rn w 
0.5, which we adopt for our calculations. The accre- 
tion mode depends on the relation between R m and the 
Keplerian co-rotation radius R co — (GM X / 1 A^v 2 ) 1 / 3, . 
Roughly speaking, for R m < R co , matter can accrete 
steadily onto the star; for R m i? co , centrifugal acceler- 
ation of infalling matter by the NS magnetosphere slows 
but does not entirely inhibit accretion, a regime we re- 
fer to as the quasi-propeller mode; and for R m > J? co , a 
true p ropeller effect as described bv llllarionov fc Sunvaevl 
(|1975|) becomes energetically possible, ejecting matter 
from the system and thereby preventing accretion. For 
AMSPs, the transition zone around R m from a Keplerian 
to a co-rotating flow is wide, so changes in the mode of 
accretion are likely to be gradual. 

Magnetohydrodynamic simulation, analytic models, 
and observational evidence all indicate that the transi- 
tion from steady accretion to a quasi-propeller regime 
occurs around R m ~ 0. 7 R rr> . T h e mag netohydro- 
dynamic simulations of ILong et "all (|2005|) identified 
this ratio as the critical magnetospheric radius below 
which v > and above which v < 0. The ana- 
lytic model of Rap paport et al.l (|2004l) produces a sim- 
ilar threshold between spin-up and spin- downQ Fi- 
nally, the well-constrained distance and magnetic field 
of SAX J1808.4— 3658 allow us to calculate its R m from 
its observed flux using eq. (f20)). The knee in the de- 
cay of that source (along with other phenomena indi- 
cating the ons et of a quasi-propeller state) occurs when 
R m = 0.6 R co (jHartman et al.ll2009b[ ). 

The magnetospheric radius implied at the knee of the 
IGR J00291+5934 outburst light curves is also consistent 



11 The ratio given in Rappaport ct al. (2004) is R m fs 1.87 R co , 
but this relation is based on a formal definition of R m that is a 
factor of 2.2 greater than ours. Adjusting for this difference gives 
R m ^ 0.85 Rco- 
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with this value. The long-term quiescent spin-down of 
this source suggests a magnetic dipole moment of /i = 
1.1 x 10 26 G cm 2 , which we derive in section [772] From 
this field strength and the critical M x derived in eq. ([15]). 
eq. (|2"U|) gives us 



10 



R 



m.kncc 



Rc 



16 km 
23.6 km 



0.7. 



(21) 



This figure is in excellent agreement with theoretical ex- 
pectations and our observations of SAX J1808.4— 3658. 
It further reinforces the validity of our assumption that 
the light curve knee marks the transition into a quasi- 
propeller accretion state. 

5.3. Halting of Accretion by the Propeller Effect? 

The transition to a true propeller state, in which the 
centrifugal acceleration is capable of ejecting matter from 
the system, occurs at R m ss 1.26 R co if the transfer 
of angular momentum is perfectly efficient but inelastic 
ijRappaport et alj 12004) ; if the transfer is a lso perfectly 
elasti c, this limit becomes R m ~ 1.13 R co (|Perna et al.1 
2006) . From eq. (|20l) , these radii respectively correspond 
to f x /f x , c « 0.12 and 0.19, giving 2.5-25 keV fluxes 
of (5-8) x 10^ n erg cm" 2 s" 1 for IGR J00291+5934. 
Within or slightly below this range, we predict another 
knee in the light curve as the decay rate again increases 
and the outburst rapidly shuts off. 

Swift XRT and XMM observations following the 2008a 
outburst require the presence of such a break. The trend 
line in Figure [7J shows the exponential decay fit to the 
PCA fluxes in the range 1 > f x / fx.c > 0.12, for which 
the NS magnetosphere should slow but not entirely halt 
accretion. Its e-folding time is 1.8 ±0.3 d, as reported in 
section [3TTI After the outburst falls below f x J f x , c ~ 0.2, 
the decay accelerates. The two last RXTE observations 
in which IGR J00291+5934 is confidently detected fall 
somewhat (2-3 a) below the extrapolated exponential de- 
cay curve. 1.4 d after the last RXTE det ection, a 3 a 
Swift XRT upper limit falls well below it (jLewis et al.l 
12010ft . Finally, a deep XMM observation 11 d after the 
outburst peak and 5 d after the last RXTE detection re- 
vealed that IGR J00291+5934 had dropped to an unab- 
sorbed2-10keVflir>Hof 1.4 x 10~ 14 erg cm" 2 s" 1 . This 
flux is lower than previously reported qu i escent levels 
by a factor of ~ 3 dJonker et al.ll2005l 120081: iTorres et all 
I2008D . although Honker et al.l (|2OO50 notes that the qui- 
escent X-ray flux can vary. IGR J00291+5934 could not 
have reached this quiescent flux so quickly without an- 
other steepening of the decay: fitting a constant decay 
rate between the last PCA detections and the XMM ob- 
servation establishes a maximum e-folding time of 0.7 d. 
The actual decay rate was likely faster. 

If the majority of optical emission during outburst is 
due to X-ray reprocessing, then we should expect a knee 
in the optical light curve if and when the mass accre- 
tion rate becomes low enough to be entirely halted by 
the propeller effect. Such a break was seen in the R 
band during the 2004 outburst when the X-ray flux was 
fx/ fx.c = 0.13 ± 0.03, in excellent agreement with our 

12 Using the spectral models of £13.21 we calculate a factor of 
1.7 to conv e rt fro m the unabsorbed 2-10 keV fluxes given by 
ILewis e t al. (2010]) to the absorbed 2.5—25 keV fluxes quoted in 
this paper. 
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Fig. 7. — A comparison of the observed fluxes and predicted 
accretion regimes during the 2008a outburst. Fluxes are from the 
PCA (dots), Swift XRT (x marks), and XMM (square ). (Swift and 
XMM measurements are taken from Lewis et al. 2010.) Horizontal 
dotted lines indicate the onset of magnetic throttling (f x / fx,c — 1 
by assumption), the approximate flux below which the propeller 
effect should halt accretion entirely (fx/fx,c ~ 0.12), and the typ- 
ical flux observed by Chandra during quiescence (see text for ref- 
erences). The solid gray line shows an exponential fit to the light 
curve during the magnetic throttling regime, and the dashed line 
extrapolates this fit to later times. Note that the XMM flux falls 
far below this trend line. There must be another break in the 
light curve around the time that the source falls below the PCA 
detection level. 



predictions (jTorres et al.ll2008h . Gaps in the optical cov- 
erage of the 2008a and 2008b outbursts around the ex- 
pecte d times of breaks p revented their detection or exclu- 
sion (jLewis et al.M2O10t F. Lewis & D. M. Russell, priv. 
comm.). 

5.4. Accretion Torques during Outburst 

Changes in the accretion state should also change the 
X-ray timing properties of IGR J00291+5934. At accre- 
tion rates above the onset of the quasi-propeller state, the 
inf ailing matter should spin up the NS. When accretion 
enters the quasi-propeller state, which we suggest corre- 
sponds to the knee in the X-ray light curve, the transfer 
of angular momentum to the NS becomes inefficient or 
negative. 

The 2004 outburst of IGR J00291+5934, which had a 
peak flux of roughly twice what was seen during 2008, 
showed a clear spin- up of 8.4(6) x 10~ 13 Hz s _1 , which 
(jBurderi et al.ll2007l) attributed to accretion torque. Un- 
der the assumption that the light curve knee approxi- 
mately corresponds to the spin-up / spin-down equilib- 
rium point, this result is expected. 

In contrast, the v measurements for the 2008 outbursts 
were poorly constrained: 10(10) x 10~ 13 Hz s _1 and 
4.5(2.5) x 10~ 13 Hz sr 1 respectively. At 1.0 er and 1.8 a 
significance, it is possible that these marginal spin-ups 
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reflect timing noise rather than any change in the spin 
of the NS. If they are real, the frequency at the end 
of the 2008a outburst would be 0.5(2) /xHz higher than 
the frequency at the beginning of the 2008b outburst, 
requiring a mean spin-down of —1.7(6) x 10 -13 Hz s _1 
during the intervening 30 d quiescence. Spin-ups dur- 
ing the 2008 outbursts would also increase the implied 
long-term spin-down. 

Yet the 2008b outburst, which never rose above / X)C , 
was likely spent entirely in the quasi-propeller state or 
at the limit of its onset. We should therefore expect 
the NS to be spun down during this outburst. This is 
at odds with its 1.8 a spin- up measurement. Given the 
low significance of this detection and the otherwise solid 
evidence that the knee at f XyC reflects a transition into 
a quasi-propeller state, it is likely that this v is due to 
timing noise rather than a change in the NS spin. The 
case of the 2008a outburst is similar: because it spent 
little time («2.5 d) above f x ,c, a large spin-up is not 
expected. 

6. MAGNETIC FIELD LIMITS FROM PULSATIONS 

The presence of accretion-powered pulsations across a 
nearly two orders of magnitude in flux constrains the 
magne tic field strength of the NS. lPsaltis fc Chakrabartyl 
(1999) used the similarly wide range of fluxes with 
detectable pulsations from SAX J1808.4— 3658 to de- 
rive limits that were compatible with the mag- 
netic field implied by the spin-down of that source 
(jHartman et al. 2009a). Here we apply their arguments 
to IGR J00291+5934. 

The detection of pulsations when IGR J00291+5934 
was at its peak flux indicates that the magnetic field 
must be strong enough to columnate the accretion flow 
above the NS surface even when the accretion rate is 
at its maximum. Setting the magnetospheric radius R m 
from eq. (|20|) equal to the NS radius R x and solving for 
the magnetic dipole moment fi gives a lower limit: 



H > 1.6 X 10'° G cm [ 



-7/4 



m x v^f Rx y/v d 



2.3 M Q J ViOkm/ \3kpc 



Cbol ' fx 



2.54 • 11.1 x 10- 10 erg cm" 2 s- 1 



1/2 



(22) 



For this conservative limit, we have taken extreme values 
of the magnetospheric constant k m , the NS mass and 
radius, and the distance to the source. 

Conversely, at the lowest fluxes with detectable pul- 
sations, the magnetic field cannot be so strong that it 
causes the ejection of matter due to the propeller effect. 
This propeller regime t urns on at around R m =1.3 R co 
(|Rappaport et alj |2004). Solving for fj,, we get an upper 
limit: 



(i < 1.2 x 10 27 G cm 



M x 



2.3 M, 



1/3 



02 

R x 



-7/4 



1/2 



(599 Hz) 



-7/6 



15 km 



5 kpc 



Cbol ' fx,: 



1/2 



2.54 • 0.66 x 10" 1U erg cm" 2 s" 



(23) 



Again, we have adopted parameter values that give a 
conservative estimate of this limit. 

7. IMPLICATIONS OF THE LONG-TERM SPIN-DOWN 

The measured 2004 and 2008 frequencies give a mean 
long-term spin-down of — (4±1) x 10 -15 Hz s _1 during the 
3.7 yr quiescence. This large v range is principally due 
to the uncertainty of whether the 2008 outburst spin-ups 
are real. If they are not, as suggested by Section l5T4l then 
the long-term spin-down is at the bottom of this range: 
-(2.8 ±0.5) x 10~ 15 Hz s _1 . For the rest of this section, 
we will use this more conservative value. Adopting the 
higher figure would only strengthen our conclusions here. 
The resulting spin-down luminosity is E = —Ait 2 Ii>i> = 



7 x 10 34 erg s 



This figure is eight times greater 



than the spin-down l uminosity of SAX J1808.4— 3658 
(jHartman et al.ll2~009a| ). due to both v and v being higher 
for IGR J00291+5934. We consider three possible mech- 
anisms for this spin-down. 

7.1. The Propeller Effect 

In section l4~3l we discussed the centrifugal expulsion of 
matter by the NS magnetosphere, commonly known as 
the propeller effect. This ejection of matter will produce 
a spin-down of the NS; however, it is unlikely that it is 
sufficient to explain the large frequency change observed 
between the 2004 and 2008 outbursts. 

The propeller effect will eject matter from near the 
magnetospheric radius, producing a torque of 



Ap r0 p s» M^GMxRm) 1 / 2 



(24) 



if matter is ejected at a rate of M cy Centrifugal ejection 
occurs if R m > 1.3 R co ~ 30 km. R m varies weakly with 
M e j, so it will not be much greater than this value. The 
Av = 0.32(4) /iHz frequency difference between the end 
of the 2004 outburst and the beginning of 2008a would 
then require a total mass ejection of ^4 x 10~ n Mq dur- 
ing the 3.7 yr quiescence. The associated rate is an order 
of magnitude greater than the expected mass transfer 
rate of3x 10~ 12 Mq yr -1 from the low-mass companion 
assuming that transfer is driven by gravitationa l radia- 
tion from the binary orbit (Galloway et al. 2005). 

7.2. Magnetic Dipole Spin-down 

Magnetic dipole spin-down is a more likely cause. Pul- 
sar magnetosphere simulations bv lSpitkov skv (2006]) give 
a magnetic dipole torque of N = — /j 2 (2nv / c) 3 (1+siia 2 a). 
Here a is the colatitude of the magnetic pole, which we 
assume to be small (~ 15°) based on the highly sinu- 
soidal pulse profil e (e.g.. iPoutanen fc Beloborodovll2006l : 
lLamb et al.ll2009T ). If this torque accounts for the en- 
tirety of the quiescent spin-down, it would require a 
magnetic dipole of (9.4 ± 0.8) x 10 25 G cm 3 . For a 
NS radius of 10 km, the corresponding surface field is 
B = 2/ii?^ 3 = 1.9 x 10 8 G at the magnetic poles. This 
field strength is consistent with the limits derived in the 
previous section from magnetospheri c arguments. It also 
agrees with the B < 3 x 10 8 G limit of lTorres et"aH (|2008[ ) 
based on the quiescent luminosity. Finally, the long-term 
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Fig. 8. — A comparison of the AMSPs (stars), Fermi pul- 
sars (black points), and radio pulsars without 7-ray emission 
(gray points). Pulsars in binary systems are circled. The 
top plot shows the standard P-P diagram. The two AM- 
SPs with known quiescent spin-downs, SAX J1808.4— 3658 and 
IGR J00291+5934, are located in the lower left corner, among 
the radio millisecond pulsars. Lines of constant ma gnetic field 
were calculated using the formula of ISpitkovskvl (I2006T) . The bot- 
tom plot shows the spin-down luminosities, normalized by dis- 
tance squared. These values are proportional to the spin-down 
flux that would be incident if all the sources r adiated with perfec t 
efficiency. Fermi pulsar parameters are from lAbdo et afl pOlOh : 
radio pulsar parameters are from the Australia Tel escope Na- 
tional Facility pulsar catalog (Manchester ct al. 2005), online at 
http://www.atnf.csiro.au/research/pulsar/psrcat Pulsars in 
globular clusters have been omitted. 



spin-down of IGR J00291+5934 places it among the ra- 
dio millisecond pulsars on the P-P diagram, suggesting 
a common spin-down mechanism. 

The high spin-down luminosity of IGR J00291+5934 
makes it a good candidate for pulsation searches dur- 
ing quiescence. The discovery of spin-powered pulsations 
from a quiescent AMSP would provide a final and elu- 
sive missing link between the AMSPs and radio millisec- 
ond pulsars. 7-ray pulsation searches may be even more 
promising. The first eight millisecond 7-ray pulsars de- 
tected by the Fermi Large Area Telescope (LAT) convert 
their spin-down luminosities into 7-ray emission at high 
efficiencies (6-100%) and emi t pulsations into wide open- 
ing angles (|Abdo et al.ll2009f ). If Fermi does not detect 
pulsations from IGR J00291+5934, it would require a 
lower dipole spin-down to 7-ray emission efficiency from 
this source, or a different spin-down mechanism entirely. 



Angular momentum loss through gravitational wave 
emission has been suggested as a way to explain the 
absence of ve ry rapid l y spinning (>730 Hz) m i llisec- 
ond pulsars (|Bildstenl 119981: IChakrabartv et al.l 120031: 
IChakrabartvl [2005) . The highly nonlinear ~i/ 5 depen- 
dence of the gravitational wave torque on spin rate v 
means that this torque would dominate at the very high- 
est spin rates but be negligible for slower spins, with 
a rather sharp transition. While there is no evidence 
that this mechanism is impor tant at spins as slow as 
400 Hz (jHartman et al.l I2008P ). IGR J00291+5934 is a 
strong candidate since it is the most rapidly spinning 
AMSP. However, as in the case of SAX J1808.4-3658 
(jHartman et al.ll2"008D , we have already found that mag- 
netic dipole torques due to the known magnetic field 
strength of the pulsar likely accounts for most of the spin- 
down in quiescence, suggesting that gravitational wave 
torques are unimportant even at 599 Hz. 

For IGR J00291+5934, the quiescent spin-down places 
an upper limit on the neutron star's mass quadrupole 
moment of^l 



Q < 1.2 x 10 



I 



10 45 g cm 2 



1/2 



(599 Hz) 



-5/2 



2.8 x 10~ 15 Hz s- 1 



1/2 



g cm 



(25) 



7.3. Gravitational Radiation Torque 



or an ellipticity of Q/I < 10 9 for moment of inertia 
/. The rapid spin of IGR J00291+5934 makes this limit 
approximately an order of magnitude more stringent 
than the similarly der ived limit for SAX J1808.4-3658 
(jHartman et al.l I2008D . It is appreciably lower than 
the upper limits on NS ellipticity from direct gravita- 
tional wave detectors: targetted LIGO searches of ra- 
dio millisecond pulsars give a strongest upper limit of 
Q/I < 7 x IP" 8 for the nearby pulsar PSR J2124-3358 
(Ab bott et"aTl[2010l ). 

The spin-down of IGR J00291+5934 begins to test pre- 
dictions of the expected quadrupole moments of accret- 
ing neutron stars. Accretion could result in a mass 
quadrupole through a variety of mechanisms. Electron 
capture due t o higher tempe r atures base of the accre- 
tion columns (Bil dstenl Il998t lUshomirskv et al.ll2000l ), 
hydrodynamic tur bulence due to accretion -induced dif- 
ferential rotation (Mclatos & Peralta 201(|), and the ac- 
cumulation of "m ountains" of accreted material (e.g., 
IHaskell et al1l2006l) all could produce ellipticities of up to 
Q/I ~ 10~ 8 ; magnetic confinement of the accreted ma- 
terial could allow an even high er limit (jMelatos &: Payne] 
120051 : iPavne fc Melatosl 120061 ). The absence of such a 
high ellipticity in IGR J00291+5934 does not rule out 
any of these mechanisms, as their prcdicitions depend on 
unknown NS parameters (e.g., the maximum strain up- 
held by the NS crust, or the ohmic diffusion time scale 
governing the settling of accreted material). Neverthe- 
less, we may be entering the regime in which it is possible 
to begin constraining some of these parameters. 

This upper limit for IGR J00291+5934 still does not 
exclude the possibility that gravitational wave torques af- 

13 The c oefficient of eq. d 2 5 H d iffers from the value derived using 
cq. (17) of Hartman et al. (2008). The coefficient in that paper is 
incorrect: it should read 1.4 X 10 36 g cm 2 , not 4.4 X 10 36 g cm 2 . 



16 



Hartman et al. 



feet the fastest pulsar spins, since the torque near the lim- 
iting «730 Hz spin frequency would be 2.7 times greater 
than in IGR J00291+5934 for the same Q. It is also 
feasible that gravitational wave emission is present only 
for a short time in these sources, acting only during and 
shortly after an outbu rst. Our limits in this paper and in 
lHartman et al.l (|2008l ) are valid only for the magnitude 
of a persistent quadrupole and do not exclude a larger 
quadrupole that dissipates substantially faster than the 
^3 yr outburst recurrence period. 



8. CONCLUSIONS 

We have presented a comprehensive analysis of the out- 
bursts of IGR J00291+5934 and a self-consistent expla- 
nation for the observed behavior. In contrast with previ- 
ous analyses of irradiated accretion disks in NS LMXBs, 
we posit that the irradiation never fully ionizes the disk of 
this source during the outbursts observed with RXTE in 
2004 and 2008. Instead, we suggest that the NS magnetic 
field and its interaction with the accretion disk plays a 
central role in shaping the light curves of these outbursts. 

The long-term spin-down of IGR J00291+5934 is most 
likely due to the magnetic dipole torque of a 2 x 10 8 G 
field at the NS surface. A magnetic field of this strength 
will impede the infall of matter when the accretion rate 
falls below 2 x 10~ 10 M© yr -1 , in excellent agreement 
with the critical flux at which knees are observed in the 
2004 and 2008a light curves. If these knees are due to 
disk-magnetosphere interaction, then it is unlikely that 
the disk was ever fully ionized: irradiation is not suf- 
ficient to ionize the entire disk at the lowest accretion 
rates observed with RXTE, and a second knee would 
likely be detectable if the transition from fully to par- 
tially ionized had occurred. Futhermore, the near-linear 
decay observed during the first week of the 2004 outburst 
is inconsistent with the exponential decay expected for a 
fully ionized disk but matches well with theoretical pre- 
dictions when partial ionization is assumed. 

Substantial evidence points to the onset of a "quasi- 
propeller" state at accretion rates below the light curve 
knee. By shutting off accretion more quickly than if the 
outburst was driven solely by irradiation, enough matter 
is left in the disk after the 2008a outburst to fuel the 



Abbott, B. P. et al. 2010, ApJ, 713, 671 

Abdo, A. A. et al. 2009, Science, 325, 848 

Abdo, A. A., Ackermann, M., Ajello, M., Atwood, W. B., 

Axelsson, M., Baldini, L., Ballet, J., Barbiellini, G., & TODO, 

X. 2010, ApJS, 187, 460 
Belloni, T., Psaltis, D., & van der Klis, M. 2002, ApJ, 572, 392 
Bildsten, L. 1998, ApJ, 501, L89 

Burderi, L., Di Salvo, T., Lavagetto, C, Menna, M. T., Papitto, 

A., Riggio, A., Iaria, R., D'Antona, F., Robba, N. R., & Stella, 

L. 2007, ApJ, 657, 961 
Burderi, L., Di Salvo, T., Menna, M. T., Riggio, A., & Papitto, 

A. 2006, ApJ, 653, L133 
Campana, S., Stella, L., Mereghetti, S., Colpi, M., Tavani, M., 

Ricci, D., Fiume, D. D., & Belloni, T. 1998, ApJ, 499, L65+ 
Chakrabarty, D. 2005, in Astronomical Society of the Pacific 

Conference Series, Vol. 328, Binary Radio Pulsars, ed. F. A. 

Rasio & I. H. Stairs, 279-289 
Chakrabarty, D., Morgan, E. H., Muno, M. P., Galloway, D. K., 

Wijnands, R., van der Klis, M., & Markwardt, C. B. 2003, 

Nature, 424, 42 



2008b outburst 30 d later. Additionally, the mass distri- 
bution left behind in the disk differs from the distribu- 
tion expected after a long period of quiescence, priming 
the disk for the slow rise and maximum at / X)C seen in 
the 2008b outburst. In the quasi-propeller regime below 
the knee, accretion torques should be small or negative, 
and the non-detection of spin derivatives during the 2008 
outbursts stands in contrast with the significant spin-up 
seen during the brighter 2004 outburst. Accretion in- 
stabilities associated with the magnetic throttling may 
explain the ~0.5 Hz QPO observed during the 2008b 
outburst. Finally, if the light curve knee is associated 
with the onset of the quasi-propeller effect, then there 
should be another knee near the RXTE PCA detection 
threshold beyond which the propeller effect can entirely 
eject infalling matter from the system, entirely halting 
accretion; indeed, Swift XRT observations and a deep 
XMM observation show that another light curve break 
must be present, causing the source to return to quiescent 
levels within days of falling below the RXTE detection 
threshold. Clearly there remains much work to be done, 
however: our invocation of magnetospheric effects has 
been largely qualitative or based on simplistic models, 
and more substantial analysis and modeling is necessary 
to gain a fuller understanding. 
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down between 2004 and 2008, although [Patrunol ()2010l ) 
claimed a higher significance for the spin-ups during out- 
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this paper. 



Chakrabarty, D., Swank, J. H., Markwardt, C. B., & Smith, E. 

2008, The Astronomer's Telegram, 1660 
de Jong, J. A., van Paradijs, J., & Augusteijn, T. 1996, A&A, 

314, 484 

de Martino, D., Matt, G., Mukai, K., Belloni, T., Bonnet-Bidaud, 
J. M., Chiappetti, L., Gansicke, B. T., Haberl, F., & Mouchet, 
M. 2001, A&A, 377, 499 

Degenaar, N., Wijnands, R., Campana, S., Galloway, D., Lcwin, 
W., Homan, J., Chakrabarty, D., Jonker, P., Cackett, E., & 
Miller, J. 2007, The Astronomer's Telegram, 1098 

Dubus, G., Lasota, J., Hameury, J., & Charles, P. 1999, MNRAS, 
303, 139 

Eckert, D., Walter, R., Kretschmar, P., Mas-Hesse, M., Palumbo, 
G. G. C, Roques, J.-P., Ubertini, P., & Winkler, C. 2004, The 
Astronomer's Telegram, 352 

Falanga, M., Bonnet-Bidaud, J. M., & Suleimanov, V. 2005a, 
A&A, 444, 561 

Falanga, M., Kuiper, L., Poutanen, J., Bonning. E. W., Hermsen, 
W., Di Salvo, T., Goldoni, P., Goldwurm, A., Shaw, S. E., & 
Stella, L. 2005b, A&A, 444, 15 



Double Outburst from ICR J00291+5934 



17 



Falanga, M. & Titarchuk, L. 2007, ApJ, 661, 1084 
Fox, D. B. & Kulkarni, S. R. 2004, The Astronomer's Telegram, 
354 

Galloway, D. 2008, in American Institute of Physics Conference 

Series, Vol. 983, 40 Years of Pulsars: Millisecond Pulsars, 

Magnetars and More, cd. C. Bassa, Z. Wang, A. Cumming, & 

V. M. Kaspi, 510-518 
Galloway, D. K. & Cumming, A. 2006, ApJ, 652, 559 
Galloway, D. K., Markwardt, C. B., Morgan, E. H., Chakrabarty, 

D., & Strohmayer, T. E. 2005, ApJ, 622, L45 
Gierlinski, M., Done, C, & Barret, D. 2002, MNRAS, 331, 141 
Gierlinski, M. & Poutanen, J. 2005, MNRAS, 359, 1261 
Gilfanov, M., Revnivtsev, M., Sunyaev, R., & Churazov, E. 1998, 

A&A, 338, L83 
Haberl, F. & Motch, C. 1995, A&A, 297, L37 
Hartman, J. M., Patruno, A., Chakrabarty, D., Kaplan, D. L., 

Markwardt, C. B., Morgan, E. H., Ray, P. S., van der Klis, M., 

& Wijnands, R. 2008, ApJ, 675, 1468 
Hartman, J. M., Patruno, A., Chakrabarty, D., Markwardt, 

C. B., Morgan, E. H., van der Klis, M., & Wijnands, R. 2009a, 

ApJ, 702, 1673 

Hartman, J. M., Watts, A. L., & Chakrabarty, D. 2009b, ApJ, 
697, 2102 

Haskell, B., Jones, D. I., & Andersson, N. 2006, MNRAS, 373, 
1423 

Heinke, C. O., Altamirano, D., Cohn, H. N., Lugger, P. M., 
Budac, S. A., Servillat, M., Linares, M., Strohmayer, T. E., 
Markwardt, C. B., Wijnands, R., Swank, J. H., Knigge, C, 
Bailyn, C, & Grindlay, J. E. 2010, ApJ, 714, 894 
Ibragimov, A. & Poutanen, J. 2009, MNRAS, 400, 492 
Illarionov, A. F. & Sunyaev, R. A. 1975, A&A, 39, 185 
Jahoda, K., Markwardt, C. B., Radeva, Y., Rots, A. H., Stark, 
M. J., Swank, J. H., Strohmayer, T. E., & Zhang, W. 2006, 
ApJS, 163, 401 

Jahoda, K., Swank, J. H., Giles, A. B., Stark, M. J., Strohmayer, 
T., Zhang, W., & Morgan, E. H. 1996, in Proc. SPIE Vol. 2808, 
p. 59-70, EUV, X-Ray, and Gamma-Ray Instrumentation for 
Astronomy VII, Oswald H. Siegmund; Mark A. Gummin; Eds., 
cd. O. H. Siegmund & M. A. Gummin, 59-70 

Jonkcr, P. G, Campana, S., Steeghs, D., Torres, M. A. P., 
Galloway, D. K., Markwardt, C. B., Chakrabarty, D., & 
Swank, J. 2005, MNRAS, 361, 511 

Jonker, P. G, Torres, M. A. P., & Steeghs, D. 2008, ApJ, 680, 615 

Juctt, A. M., Galloway, D. K., & Chakrabarty, D. 2003, ApJ, 
587, 754 

King, A. R. & Ritter, H. 1998, MNRAS, 293, L42 (KR98) 
Lamb, F. K., Boutloukos, S., Van Wassenhove, S., Chamberlain, 

R. T., Lo, K. H., Clare, A., Yu, W., & Miller, M. C. 2009, ApJ, 

706, 417 

Lewis, F. ., Linares, M., Russell, D. M., Wijnands, R., & Roche, 
P. 2008, The Astronomer's Telegram, 1726 

Lewis, F., Russell, D. M., Jonker, P. G., Linares, M., Tudose, V., 
Roche, P., Clark, J. S., Torres, M. A. P., Maitra, D., Bassa, 
C. G., Steeghs, D., Patruno, A., Migliari, S., Wijnands, R., 
Nelemans, G., Kewley, L. J., Stroud, V. E., Modjaz, M., 
Bloom, J. S., Blake, C. H., & Starr, D. 2010, A&A, in press 
(arXiv: 1005.1178) 

Linares, M., van der Klis, M., & Wijnands, R. 2007, ApJ, 660, 595 

Linares, M., Wijnands, R., van der Klis, M., Krimm, H., 
Markwardt, C. B., & Chakrabarty, D. 2008, ApJ, 677, 515 



Long, M., Romanova, M. M., & Lovelace, R. V. E. 2005, ApJ, 
634, 1214 

Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, 
AJ, 129, 1993 

Markwardt, C. B., Galloway, D. K., Chakrabarty, D., Morgan, 
E. H., & Strohmayer, T. E. 2004a, The Astronomer's Telegram, 
360 

Markwardt, C. B., Pereira, D., & Swank, J. H. 2007, The 

Astronomer's Telegram, 1051 
Markwardt, C. B., Swank, J. H., & Strohmayer, T. E. 2004b, The 

Astronomer's Telegram, 353 
Mclatos, A. & Payne, D. J. B. 2005, ApJ, 623, 1044 
Melatos, A. & Peralta, C. 2010, ApJ, 709, 77 
Muno, M. P., Ozel, F., & Chakrabarty, D. 2002, ApJ, 581, 550 
Paizis, A., Nowak, M. A., Wilms, J., J-L. Courvoisier, T., 

Ebisawa, K., Rodriguez, J., & Ubertini, P. 2005, A&A, 444, 357 
Papitto, A., Di Salvo, T., Burderi, L., Menna, M. T., Lavagetto, 

G., & Riggio, A. 2007, MNRAS, 375, 971 
Papitto, A., Riggio, A., Burderi, L., Di Salvo, T., D'Ai, A., & 

Iaria, R. 2010, A&A, submitted (arXiv: 1006.1303) 
Patruno, A. 2010, ApJ, 722, 909 

Patruno, A., Watts, A. L., Klein- Wolt, M., Wijnands, R., & van 

der Klis, M. 2009a, ApJ, 707, 1296 
Patruno, A., Wijnands, R., & van der Klis, M. 2009b, ApJ, 698, 

L60 

Payne, D. J. B. & Melatos, A. 2006, ApJ, 641, 471 
Pcrna, R., Bozzo, E., & Stella, L. 2006, ApJ, 639, 363 
Poutanen, J. & Beloborodov, A. M. 2006, MNRAS, 373, 836 
Powell, C, Haswcll, C, & Falanga, M. 2007, MNRAS, 374, 466 
(P07) 

Pringle, J. E. 1996, MNRAS, 281, 357 
Psaltis, D. & Chakrabarty, D. 1999, ApJ, 521, 332 
Rappaport, S. A., Fregeau, J. M., & Spruit, H. 2004, ApJ, 606, 
436 

Rcmillard, R. 2004, The Astronomer's Telegram, 357 
Rupen, M. P., Dhawan, V., & Mioduszewski, A. J. 2004, The 

Astronomer's Telegram, 364 
Shahbaz, T., Charles, P. A., & King, A. R. 1998, MNRAS, 301, 

382 

Shakura, N. I. & Syunyaev, R. A. 1973, A&A, 24, 337 

Spitkovsky, A. 2006, ApJ, 648, L51 

Spruit, H. C. & Taam, R. E. 1993, ApJ, 402, 593 

Torres, M. A. P., Jonker, P. G, Steeghs, D., Roelofs, G. H. A., 

Bloom, J. S., Casares, J., Falco, E. E., Garcia, M. R., Marsh, 

T. R., Mcndez, M., Miller, J. M., Nelemans, G., & 

Rodn'guez-Gil, P. 2008, ApJ, 672, 1079 
Ushomirsky, G, Cutler, C, & Bildstcn, L. 2000, MNRAS, 319, 

902 

van der Klis, M. 1995, in The Lives of the Neutron Stars, ed. 
M. A. Alpar, U. Kiziloglu, & J. van Paradijs, 301-330 

Wood, K. S., Ray, P. S., Bandyopadhyay, R. M., Wolff, M. T., 
Fritz, G, Hertz, P., Kowalski, M. P., Lovellette, M. N., Yentis, 
D., Bloom, E. D., Giebels, B., Godfrey, G, Reilly, K., Saz 
Parkinson, P., Shabad, G., & Scargle, J. 2000, ApJ, 544, L45 

Wood, K. S., Titarchuk, L., Ray, P. S., Wolff, M. T., Lovellette, 
M. N., & Bandyopadhyay, R. M. 2001, ApJ, 563, 246 



